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Abstract 

In the framework of the Schwinger-Dyson equation and the Bethe-Salpeter 
equation in the improved ladder approximation, we calculate the — vr^ mass 
difference on the same footing as the pion decay constant and the QCD S 
parameter (or Lio) through the difference between the vector current correlator 
Hvv and the axial- vector current correlator ILaa- We find that all the results 
can be fit to the experimental values for rather large Aqcd ~ TOOMeV which 
refiects the "scale ambiguity" . By fitting to the calculated data using the pole 
saturated form of Hvv — ^AA, we also derive masses and decay constants of 
p meson and ai meson, which we found are consistent with the experiments 
rather insensitively to the "scale ambiguity" . 
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1 Introduction 



The 7r~^ - tt^ mass difference Am^ = — m^o is an interesting quantity to measure 
an explicit breaking of tlie chiral symmetry by the gauge couphng (f/(l)em) in the 
spontaneously broken phase of the chiral symmetry. It in fact has been a prototype 
of the mass calculation of pseudo Nambu-Goldstone (NG) bosons in strong coupling 
gauge theories such as those in the technicolor theories PP and more recently in the 
little Higgs models . The sign as well as the absolute value of Am^ is an important 
issue ("vacuum alignment problem"), since the negative sign would imply that the 
vacuum would align so as to break spontaneously the gauge symmetry (f/(l)em in 
the case of pion) which explicitly breaks the spontaneously broken symmetry. Such 
a situation never happens in the real-life QCD but may do in other theories. Hence 
this type of calculation plays a central role of the model buildings. 

The first successful calculation of Am^ j3] was done by the current algebra in 
conjunction with the Weinberg spectral function sum rules saturated by the n, 
p and Oi meson poles. Somewhat more elaborate calculation on this line was done 
combined with the QCD information Recently, an effective field theory calculation 
without Oi meson has been successfully done jH], based on the Hidden Local Symmetry 
(HLS) model [3 El El at loop level jTUI HH and the Wilsonian matching [H] of the 
HLS model with the QCD (for reviews of the HLS approach, see ^lE])- Ref. jU] 
demonstrated that the HLS model is a little Higgs model with two sites and two links 
(open moose) whose locality of theory space forbids quadratic divergence in Ami 
even without ai meson. 

On the other hand, direct QCD calculation of the Am^ has never been done ex- 
cept for the lattice simulation ^H]- Although the lattice QCD is a powerful method, 
it is still important to establish alternative QCD-based methods which are more in- 
tuitive and less dependent on the power of computers. Such methods will be crucial 
to the subjects which are difficult to be studied in the framework of the lattice simu- 
lation. Most relevant subjects are the dynamical origin of the electroweak symmetry 
breaking, QCD at finite density, etc. 

In this paper, we calculate in the chiral symmetric limit the Ami simultaneously 
with the the pion decay constant /^r and the QCD S parameter (TH] (or the Gasser- 
Leutwyler parameter Liq ^7j) within the framework of the the Schwinger-Dyson 
(SD) equation and the Bethe-Salpeter (BS) equation in the the improved ladder 
approximation. We note that all these quantities are given by "generalized" Weinberg 
spectral function sum rules for the difference between the vector current correlator 
Hvv and the axial- vector current correlator Haa- fl is given by the Weinberg first sum 
rule, the S parameter by the Das-Mathur-Okubo (DM0) sum rule jlH] or the "zeroth 
Weinberg sum rule", and Am^ by the Das-Guralnik-Mathur-Low- Young (DGMLY) 
sum rule [3j or the "third Weinberg sum rule" . We then rewrite them in terms of the 
current correlators instead of the spectral functions and directly calculate Am^ on 
the same footing as /tt and the QCD S parameter through liyv — ^aa in the space- 
like momentum region. Such correlators are written in terms of the BS amplitudes 
which are given by solving the inhomogeneous BS (IBS) equation and the SD equation 
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with the improved ladder approximation in the Landau gauge both in the space-like 
region. In contrast to our case, solving the bound state spectra by the homogenous 
BS equation in the time-like region is usually difficult, since it would need analytic 
continuation of the running coupling which is not an analytic function. ^ 

This kind of method has been extensively used to investigate masses and decay 
constants of low lying mesons [201 1^ 122] in good agreement with the experiments. 
On the other hand, the same method failed to reproduce the QCD S parameter 
consistently with the experiment |23j . 

However, we find that the QCD S parameter is rather sensitive to the infrared 
(IR) cutoff parameter of the QCD running coupling. We actually reproduce the 
QCD 5* parameter consistently with the experiment by setting the IR cutoff param- 
eter in the region never investigated in the previous works. Our parameter choice 
corresponds to taking the QCD scale parameter as Aqcd — 724 MeV, while it was 
chosen as Aqcd — 500 MeV in the previous results of the decay constants and the 
masses [201 121] and of the QCD S parameter both values being considerably 
higher than the conventional Aq^j^ = 300-450 MeV obtained in the MS scheme j24j . 
This actually corresponds to exploiting the "scale ambiguity" [25j of Aqcd in the SD 
and BS equations recently emphasized by Hashimoto and Tanabashi |26j . 

We then show that 7r+ — 7r° mass difference Am^ is calculated in rough agreement 
with the experimental value for the same parameter region as that reproducing the 
QCD S parameter and /^r. Actually, although infrared dynamics is quite important 
for determination of /^r and S, it is not so important for determination of Am^. In 
other words, and S are sensitive to the infrared structure of the running coupling, 
while Am^ is not. 

We also derive the masses and decay constants of p meson and ai meson by fitting 
them to the calculated Hyv ~ ^aa within our method using the pole saturated form 
of Hyy — H^^. These quantities are also compared with experiments, which is another 
check of the validity of the calculation in the present analysis. 

The calculation presented in this paper turns out to be the ffist example which 
derives Am^ directly from QCD without depending on the lattice simulation. Once 
we have checked the reliability of such a method in the ordinary QCD (with three light 
flavors), we may apply it to the mass of pseudo NG bosons and vacuum alignment 
problem in other strong coupling gauge theories such as the large Nf QCD, walking 
technicolor, little Higgs, hot/dense QCD, etc.. 

This paper is organized as follows. In section |2] we introduce spectral function 
sum rules: The DM0 sum rule or the zeroth Weinberg sum rule for S parameter, 
the first Weinberg sum rule for the second Weinberg sum rule, and the DGMLY 
sum rule or the third Weinberg sum rule for Am^. We then rewrite them in terms of 
the current correlators Hvv — ^aa- In section El we show how the current correlators 
are obtained from the BS amplitude which is calculated from the IBS equation in the 
space-like region. In section|3]we introduce the IBS equation and show how to solve it 
numerically. Section |3] is the main part of this paper: We ffist show the result of the 

^See an exception in Ref. [T^. 
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IBS equation for QCD, and discuss the dependence of f-^, S and Am^ on the infrared 
structure of the running couphng. Then we compare them with the experimental 
values. We also show results for masses and decay constants of p and ai mesons. In 
section El we give summary and discussions on applications to the large Nf QCD and 
the electroweak symmetry breaking. 



2 Spectral function sum rules 

In this section we introduce the spectral function sum rules, which express the pion 
decay constant the QCD S parameter and Am^ in terms of Hyv ~ ^aa in the 
space-like momentum region. Here we consider only the chiral symmetric limit of 
massless three flavors [Nf = 3), the corrections of the finite quark mass being expected 
to be small for the quantities we consider. 

Let us begin with introducing the vector and axial-vector currents as 

V;{x)=ijix)T^^,i;{x), (2.1) 

A^ix) = V^(x)T"7^75V^(x) , (2.2) 

where T'^ is the generator of SU{Nf ) normalized as tr(T"T^) = and we consider 
a,b = 1, 2, 3. In the chiral symmetric limit the current correlator IIjj is given as 



'^"M ^ - j njj(g') = ^Jd'xe''^^{0\TJ;{x)JtmO), (2.3) 

iJ;ix) = V^ix),A;ix)). 



q2 



The Umezawa-Kamefuchi-Kallen-Lehmann spectral representation for the current 
correlators are expressed as 

d^x e*«" (0|rj"(x)J,^(0)|0) 



= -S^' [ ^^7373^ { {3,^ - ^) P^^^'^ - Wfis)] ■ (2.4) 

Here, we decomposed the spectral function into the spin-one part p^j\s) and spin- zero 
part pj'*(s). Since the axial- vector current couples to the massless NG boson, the 
pion, while no massless particles couple to the vector current V^, we have 

p^is) = 0, (2.5) 
pf{s) = fJis), (2.6) 

where /tt is the decay constant of the NG boson vr defined by 

(0|A;:(0)|vr^(g))=^g,/X^ (2.7) 
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Using the spectral functions, we can write down the following "generalized" Wein- 
berg sum rules: 
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where, 5* is the QCD 5* parameter (1^, which in the case of Nj 
the Gasser-Leutwyler parameter Lio [Hj as S* = — IGvrLio, (yem{= 
coupling constant of the electromagnetic interaction, and Ami the tt 
difference defined by Am^ = — m^o • Equations (|2.9|) and (j2.10p are the first and 
second Weinberg sum rules ji], respectively, and Eqs. ()2.8p is the DM0 sum rule [TH] 
or often called the "zeroth Weinberg sum rule" . Equation 1)2.111) was derived by Das 
et al. jS] and may be called the DGMLY sum rule or the "third Weinberg sum rule" . 

Of course these relations in Eqs. ()2.8|1 . ()2.9|) . ()2.10|) . and ()2.11|) make sense only 
if the integrals converge. For example, convergence of the second sum rule ()2.10j) 
requires that the difference between Hvv and Haa must satisfy 



(Q2^0o) 



0, 



(2.12) 



where is related to the space-like momentum as = —q^{> 0) . 

The asymptotic behavior of the current correlators can be calculated by the oper- 
ator product expansion (OPE) technique. In the chiral limit, the form of liyviQ'^) ~ 
^aa{Q'^) in the ultraviolet region is estimated as follows |77ll^ : 



nvy(Q')-n^A(g' 



.oo) 47r(iV2-l) as{qqf 



(2.13) 



up to logarithm. This means that UyviQ'^) ~ ^aa{Q^) indeed satisfies the condition 
()2.12|) and the relation in Eq. ()2.10|) makes sense. When the integral in the second 
sum rule converges, those in the first and zeroth sum rules are also convergent. 

Now we rewrite the zeroth and the first sum rules in terms of the current correlators 
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From Eq. ()2.16p . we can see that the asymptotic behavior ()2.13|) also guarantees con- 
vergence of the third sum rule. 

In Ref. jSSj the pion decay constant /^r and the QCD S parameter were calculated 
by the use of the above formulas. However, Am^ was not calculated so far. So the 
calculation presented in this paper turns out to be the first example which derives 
Ami directly from QCD without depending on the lattice simulation. Moreover, 
there is another interest regarding Am^ related to the structure of the QCD vacuum, 
as we stressed in the Introduction. When we switch off the electromagnetic interaction 
and set masses of quarks to be zero, pions are identified with the exact NG bosons 
associated with the spontaneous breaking of SU{2)l(^ SU(2)ii chiral symmetry down 
to SU{2)v symmetry. The existence of the electromagnetic interaction explicitly 
breaks the SU{2)l® SU{2)fi chiral symmetry, which makes vr"*" and vr" be the pseudo 
NG bosons. On the other hand, 7r° remains massless since the photon does not 
interact with 7r°. The interesting point here is whether Am^ = — m^o (= "02^+ 
in the chiral limit) becomes positive or negative, which is called "vacuum alignment 
problem". Negative Am^ means that fluctuation of tt"*" field around (vr^) = is 
unstable and the vacuum with {tt'^) = is not a true vacuum. If this is the case, vr"*" 
has non-zero vacuum expectation value and f/(l)cm symmetry is broken. We know 
that, in the real world, Am^ is positive and the vacuum with (7r+) = is the true 
vacuum. U{l)em symmetry is not broken in the real world. This is quite nontrivial 
fact resulting from the nonperturbative dynamics of the strong interaction. So it is 
interesting to investigate whether we can reproduce positive Am^ which is realized 
in the real-life QCD. 

3 Current correlators from BS amplitudes 

In the previous section we have written down the QCD S parameter, /^r, and Am^ 
in terms of the current correlators. Then, in this section, we show how the current 
correlators are obtained from the BS amplitude which will be calculated from the IBS 
equation. 

To derive properties of hadrons as boundstates, it is straightforward to perform 
calculations directly in the time-like momentum region. However, it is difficult to 
solve the BS equation and the SD equation in the time-like region since we have to 
carry out the analytic continuation of the running coupling from the space-like region 
to the time-like region. In the case of QCD, the one-loop running coupling is not an 
analytic function, so that we have to approximate it by some analytic function in order 
to perform the calculations in the time-like region. However, here we need the BS 
amplitude only for the space-like region in order to calculate the current correlators 
for Am^, and the QCD S parameter. 

The BS amplitude x^"^^ {J = ^) ^) is defined in terms of the three-point vertex 
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function as follows: 



A" 



/' 



-ipr. 



xi'^{p;q,e) = e" j dW«-(0|T V^«,;(r/2) (-r/2) |0), 

(3.1) 

where is the total momentum of the fermion and anti-fermion and p^ is the relative 
one. is the polarization vector defined by e-q — 0, e-e = —1, and (/, /'), (i, j), (a, (3) 
are flavor, color and spinor indices, respectively. Closing the fermion legs of the above 
three-point vertex function and taking the limit r ^ 0, we can express the current 
correlator in terms of the BS amplitude as follows: 



njj(g' 



34" 



d^p N, 



tr 



i(27r)4 2 
where 

Gir^ = 7m, G\f^ = 7.75, (3.3) 

and Ai'c = 3 is the number of colors. In the above expression wc averaged over the 
polarizations so that IT jj(g^) does not depend on the polarization. 



{e.G^'^)x^'\p-q,e) 



(3.2) 



We expand the BS amphtude XafiiP'i ^) terms of the bispinor bases Vl"^' and 
the invariant amplitudes Xi^^ ^ 



(J) 



x^'^^ (p; q, e)] = E [r^^ (p; q^ ^)] xf^ (p; q) , 



1=1 



(3.4) 



where = g^^/ The bispinor bases can be chosen in a way that they have the 
same properties of spin, parity and charge conjugation as the corresponding current 
J^{x) has. We adopt the following bispinor bases for the vector vertex: 

rr = ^ rr) = i[^^](p.g), r^^li^J], rf^ = ^[/,^,^] (3.5) 

where [a, 6, c] = a[6, c] + &[c, a] + c[a, 6]. For the axial- vector vertex we use 



{A) 



/^75, r. 



{A) 



2 -^i^mi.. ri-^^ = 2[^^]b-^)75, 



3! 



= ii^-p) (p-q) 75, r^"^^ = ^[/^, • p) (p • 75- 



<A) 



(3.6) 



Prom the above choice of the bases, we can easily show that all the invariant ampli- 
tudes Xj-"^^ are the even functions of (p • q) using the charge conjugation property of 
the current. 

In the present analysis we fix the frame of reference in such a way that only the 
zero component of the total momentum q^ becomes non-zero. Furthermore, we study 
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the case where q'^ is in the space-hke region. Then, it is convenient to parameterize 
the total momentum g'^ as 

g'^ = (zQ, 0,0,0). (3.7) 

For the relative momentum p^^, we perform the Wick rotation, and parameterize it 
by the real variables u and x as 



p ■ q = —Q u , p^ = —u^ — 



Consequently, the invariant amplitudes xl^"* become functions in u and x: 



{3.i 



(3.9) 



From the charge conjugation properties for the BS amplitude x^"^^ and the bispinor 
bases defined above, the invariant amplitudes Xi"'^(w,a;) are shown to satisfy 



Using this property of the invariant amplitudes, we rewrite Eq. ()3.2j) as 
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xj^^ (u, X] Q) - y Xe"^^ {u, x; Q) 



(3.10) 



(3.11) 



(3.12) 



Here, we used the expanded form of the BS amplitude shown in Eq. ()3.4|) and carried 
out the three dimensional angle integration. 

From Eqs. (j3.11|) and (j3.12p . the quantity liyv ~ ^aa is expressed as 
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+ j{x^^\u,x;Q) + xi''\u,x;Q) 



. (3.13) 



We note that, although either Hvv or Haa is logarithmically divergent quantity, the 
difference Hvv — ^aa becomes finite due to the cancellation of the divergence ensured 
by the chiral symmetry. 



4 Inhomogeneous Bethe-Salpeter equation 

In this section we introduce the inhomogeneous Bethe-Salpeter (IBS) equation from 
which we calculate the BS amplitude defined in the previous section. We also show 
the numerical method for solving the IBS equation. 
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Figure 1: A graphical expression of the IBS equation in the (improved) ladder approx- 
imation. 

4.1 IBS equation 

The IBS equation is the self-consistent equation for the BS amphtude x'''^\ ^ind it is 
expressed as (see Fig. [T] for graphical expression) 

T{p-q)x^'\v;q.e) = e ■ G^'^ + K{p;k) * x^''\k;q,e). (4.1) 

The kinetic part T is given by 

T{p; q) = Sp\p + q/2) ® S^\p - q/2) , (4.2) 

where Sp is the full fermion propagator iSp^{p) = — S(p). (Note that wave function 
renormalization factor A{p) becomes unity when we use the Landau gauge.) The BS 
kernel K in the improved ladder approximation is expressed as 

where g{p, k) is the running coupling of QCD whose explicit form will be shown later. 
In the above expressions we used the tensor product notation 

{A(»B)x = AxB , (4.4) 

and the inner product notation 

K{p; k) * x^'\k- q,e)=J K{p, k) x{k- q) . (4.5) 

The mass function of the quark propagator is obtained from the SD equation: 

j:{p) = K{p,k)*tSF{p). (4.6) 

It should be stressed that we must use the same kernel K{p, k) as that used in the 
IBS equation for consistency with the chiral symmetry [HTH IHT| 1^ \6'6\. 
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4.2 Numerical method for solving the IBS equation 

In this subsection we transform the IBS equation in Eq. ()4.H) into the form by which 
we can solve it numerically. 

First, we introduce the conjugate bispinor bases defined by 

rrHp;?,e)=7orrHp*;g,e)S. (4.7) 

Multiplying these conjugate bispinor bases from left, taking the trace of spinor indices 
and summing over the polarizations, we rewrite Eq. 1)4.11) into the following form: 

oTT J —oo Jo 

(4.8) 

where the summation over the index j is understood, and 

4'^ = Etr[rS"^(p;g,e)(e-G(-^))] , (4.9) 

e 

tI'\u,x) = Y.^T[Ti'\p;q,e)T{p;q)T'j'\p;q,e)], (4.10) 

e 

Klpiu,x;v,y) = j'^dcosO J2^T[T\'\p;q,e)Kip,k)Tp\k;q,e)] , (4.11) 

with the real variables v and y introduced as 

k-q = —vQ, k ■ p = —uv — xy cos9. (4-12) 

Here 6 is the angle between the spatial components of and fc^. 

Using the property of Xi^^ in Eq. (jH.lOj) . we restrict the integration range as f > 0: 



dvKij{u,x;v,y)Xj \v,y) = / dv [Kij{u, x]v,y) + Kij{u, x; -v,y)]xj iv,y). 

(4.13) 

Then, in the following, we treat all the variables u, x, v and y as positive values. 

To discretize the variables u, x, v and y we introduce new variables U, X, V and 
Y as 

TT V 

u = e , X = e , 
= , y = e^ , (4.14) 

and set ultraviolet (UV) and infrared (IR) cutoffs as 

U,V e [Xu,Au], X,Y e [Ax, Ax]. (4.15) 

We discretize the variables U and V into Nbs,u points evenly, and X and Y into 
Nbs,x points. Then, the original variables are labeled as 

U[i^] = exp [\u + Dulu] , = exp [Ax + Dxlx] , 

V[j^] = exp [Ac/ + Duly] , y[iy] = exp [Ax + Dxly] , 
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where lu, /y = 0, 1, 2, ■ ■ ■ {Nbs,u — 1) and Jx, /y = 0, 1, 2, ■ ■ ■ {Nbs,x — !)• The mea- 
sures Du and Dx are defined as 



D 



Au-Xu 



u 



D 



Ax- A; 



X 



Nbs,u — 1 Nbs,x — 1 

As a result, the integration is converted into the summation: 



/ dy dv 

Jv>0 



DuDy J2 V 
IvJy 



(4.16) 



(4.17) 



In order to avoid integrable singularities in the kernel K{u, x; v, y) at {u, x) = {v, y), 
we adopt the following four-splitting prescription 



Kij{u,x,v,y) 



1 



Kij{u,x,v+,y+) + Kij{u,x,v+,y^) 
+ Kij{u, X, v-,y+) + Kij{u, X, V-,y-) ], (4.18) 



where 



v± = exp 



V± 



D 



u 



y± = exp 



Y± 



D 



X 



(4.19) 



Now that all the variables have become discrete and the original integral equation 
(14. ip has turned into a linear algebraic one, we are able to deal it numerically. 



5 Results and discussion 

In this section, we show the results of calculations and give some discussions. First, 
we discuss the dependence of the results on the infrared structure of the running 
coupling. Second, we compare the results with experimental values. 

For solving the IBS and SD equations, we have to fix the form of the running cou- 
pling g'^{p, k) which appears in the IBS and the SD equations. We use the solution of 
the renormalization group equation for the QCD running coupling with one-loop ap- 
proximation. We regularize the infrared divergence of the one-loop running coupling 
by introducing the IR cutoff parameter tp (> 0) as follows: 

a{p,k)^^^^ = ao j-— t = \n[ipl + kl)/Alcn\ (5.1) 

where Oq = 127r/(llA'"c — 2A^j) with Nj being the number of flavors. As for the 
argument of the running coupling g'^{p, k), we used the angle averaged form p|; + k'^ 
so that we can analytically carry out the angle integration in the SD and the IBS 
equations. 

In Fig. El we plot the solutions of the SD equation for several values of tp- For 
smaller value of tp, the value of mass functions S(a;) in the infrared energy region 
becomes larger. This is natural because smaller value of tp means larger running 
coupling in the infrared region. 
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Figure 2: Mass functions obtained from the SD equation for several values of tp. 

When we solve the IBS equation, we use the following parameters: 

[At;,Ac7] = [-5.5,2.5], (5.2) 

[Ax, Ax] = [-2.5,2.5], (5.3) 
Nbs,u = Nbs,x = 30. (5.4) 

These parameters are chosen so that the dominant supports always lie within the 
energy region between UV and IR cutoffs for all values of tp which we use in the 
present analysis (tp = 0.04-0.5). 



5.1 tp dependence of Hyv — ^aa 

In Fig. El we plot the resultant values of UyviQ'^) — ^aa^Q^) for several values of 
tp. Horizontal axis shows the square of the Euclidean momentum Q. In this figure, 
HvviQ"^) ~^aa{Q'^) and are normalized by Aq^^, and plotted as the dimensionless 
quantities. As we mentioned before, either UyviQ'^) or IIaa{Q'^) is logarithmically 
divergent. These divergences are expected to be canceled in UwiQ"^) — ^aa^Q"^) 
ensured by the chiral symmetry. Figure El shows this cancellation actually occurs 
in the present calculation. ^ From the results shown in Fig. El we can see that 
UvviQ^) ~ ^aa{Q^) is dependent on the value of tp. In other words, UwiQ^) ~ 
^aaIQ"^) is sensitive to the infrared structure of the QCD running coupling. 

Once we obtain UvviQ"^) — ^aa^Q"^), we can calculate the QCD S parameter, 
and Ami by using Eqs. ()2.14p . ()2.15j) and ()2.16|) . In table HI we list the resultant 

^If the quadratic divergences were not canceled, Tlyv — ^aa would become logA^ ~ 

O.lAg^^. Figure 13 shows that Uyy — Haa is about O.OSAg^^^ at biggest, so that the divergences 
are actually canceled with each other. 
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Figure 3: Resultant values of UvviQ^) — IIaa{Q^) for several values oitp- Horizontal 
axis shows the square of the Euclidean momentum normalized by the scale Ag^^,. 
^vv{Q^) — Haa{Q^) in the vertical axis is also normalized by Agj^^, 
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0.04 


0.05 


0.1 
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0.3 


0.5 
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0.33 


0.33 


0.35 


0.39 


0.43 


0.47 


fn /J^QCD 


0.128 


0.132 


0.145 


0.162 


0.171 


0.171 


^< I^QCD 


0.00201 


0.00203 


0.00212 


0.00215 


0.00212 


0.00201 



Table 1: Resultant values of the QCD S parameter, and Am^ for several values of 
tp. DimensionfuU quantities are normalized by Aqcd and expressed as dimensionless 
quantities. 



values of these quantities for several values of tp- Here, we normalized dimensionfull 
quantities by Aqcd, and expressed them as dimensionless quantities. These results 
show that the QCD S parameter and /tt/Aqcd o.re sensitive to the value of IR cutoff 
parameter tp while ^rn^/ Kq(jj~, is less sensitive to it. This is natural because the QCD 
S parameter and /^r are directly related to the infrared quantities as can be easily 
seen from Eqs. fl2.14|l and ()2.15j) . while the infrared dependence of the integration of 
nFy(Q^) - HaaIQ^) in Eq. (E:T!)|1 is compensated by fl{= nyy(O) - Haa^^)) in the 
denominator. 

Tabled shows that, when we take the value oitp smaller, the value of fn/AqcD 
becomes smaller. Here, smaller value of ti? means larger value of infrared running 
coupling, or, in other words, strong infrared dynamics. At first sight, this behavior 
seems strange because is the order parameter of the chiral symmetry breaking, and 
it is expected to reflect the magnitude of the chiral symmetry breaking. One might 
think that /tt/Aqcd becomes larger for smaller value of tp. However, it is not the 
case. In the present analysis, /,r is not necessarily proportional to the magnitude of 
the chiral symmetry breaking because QCD with Nf = 3 in the vacuum is very far 
from the chiral phase transition point. To understand this behavior, let us see the 
Pagels-Stokar (PS) formula 



3. 



As shown in Fig. |2l the mass functions S(x) are almost constant for < x = < 
Aqcd suddenly drops at x = I^qcd all values oi tp. So we can cut the 
integration in Eq. ()5.5j) at x = I^qcd and drop the derivative of mass function. We 
also drop x in the denominator since x satisfies x -C S(x) for the relevant integral 
region (0 < x < Aq^^). Then we can approximate Eq. ()5.5|) as 

Nr. r^QCD ^ 1 K%CD 



^ 4^2 X '^'^'^£2(0) S2(0)- ^^-^^ 

From this, we can easily understand that, when infrared dynamics becomes strong, 
i.e., S(0) becomes large, becomes small. 

^We checked that the difference between calculated from the IBS equation and the PS formula 
is less than 10 % for all values oitp. We also checked that qualitative feature oi tp dependence of 
calculated from the PS formula is not different from /^r obtained from the IBS equation. 
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tp 


0.04 


0.05 


0.1 


0.2 


0.3 


0.5 


Exp. 


Ami (MeV2) 


1050 


1003 


855 


698 


620 


585 


1261.1 


S 


0.33 


0.33 


0.35 


0.39 


0.43 


0.47 


0.32 ± 0.04 


U (MeV) (input) 


92.4 


92.4 


92.4 


92.4 


92.4 


92.4 


92.4 



Table 2: Resultant values of Am^, S and /^r for several values of tp- Here, we used 
/tt = 92.4 MeV as an input to introduce the physical energy scale. Experimental values 
of Ato^ and f.^ are given in Ref. and S is obtained from the value of the Liq 

given in Ref. [Tlj through S = -~16n Lio{^) 



1 



In^ 



For comparison, we consider the decay constant near the chiral phase transition 
point. In this case, S(a;) becomes very small and Eq. ()5.5p is well approximated by 

From this, we can see that if the running coupling, and then the mass function be- 
comes smaller, f^^ also becomes smaller. Thus near the critical point, order parameter 
/,r is actually proportional to the magnitude of the chiral symmetry breaking. We can 
actually see this behavior in the large Nf QCD near the critical point (see Figs. 4, 5 
in Ref. 



5.2 Comparison with experiments 

In this subsection, we compare our results with experiments. So far in this paper, 
all the dimensionfull quantities were normalized by Aqcd and dealt as dimensionless 
quantities. Here, we introduce the "physical energy scale" by setting Aqcd in a 
way that it reproduces the experimental value of Since f-,, is sensitive to the 
infrared parameter, fixing implies that Aqcd in turn appears sensitive to the 
infrared parameter and so does the Ami which was shown insensitive to the infrared 
parameter as far as Aqcd is fixed. We use = 92.4 MeV as an input to fix Aqcd in 
the present analysis.^ 

Table|21shows comparison of the resultant values of Am^ and S for several values of 
tp with experimental values. Am^, S and /tt are directly calculated from YlwiQ"^) — 
IlAAiQ"^) by using Eqs. (jTTH|l . ^IM and (jTTH|l . From table we can see that if 
we change the value of tp, values of Am^ and S move to opposite directions: When 
we decrease the value of tp, Ami becomes large, while S becomes small. For the 
choice tp = 0.5, which was adopted in Ref. [SEl 1201 as a good regularization in 
the calculations of the masses and decay constants of p and ai mesons, we obtained 
Ami ~ 585 MeV^ and S = 0.47. Both of these values are far from their experimental 
values Ami = 1261.1 MeV^ and S = 0.29 - 0.36. 

^In the chiral symmetry limit of three flavors to„ = = rris — 0, the pion decay constant is 
estimated as /^r = 86.4 ± 0.26 MeV ^1). In this paper we use instead the physical value only for 
fixing Aqcd as a reference scale, which is up to the "scale ambiguity" in the SD and BS approach. 
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However, as we decrease the value of tp, both Am^ and S move toward their 
experimental values. At tp = 0.04, they become 



Ami = 1050 Me V^ (5.8) 
S = 0.33, (5.9) 

in good agreement with the experimental value of S and in rough agreement with 
that of Am^. In the case of tp = 0.04 we take 

AQc;D = 724MeV (5.10) 

in order to set = 92.4 MeV. This value is much larger than the conventional value 
of Aqcd in the MS scheme ^ : 

Ag)^j3 = 300 - 450 MeV . (5.11) 

Such a problem already occurred in the previous works jSHl 1201 121] , 

Aqcd - 500 MeV, (5.12) 

which roughly corresponds to tp = 0.5 in the present analysis. This reflects the "scale 
ambiguity" as was recently pointed out in Ref. [26j where an interesting method 
was proposed to solve it by using the effective couphng. 

From these results, we conclude that, when we calculate f^, Am^ and S by using 
the BS equation with the improved ladder approximation, the infrared cutoff parame- 
ter tp should be taken smaller than values used so far in the previous works ^Hl 12011211 • 
This means that the running coupling in the infrared energy region should be taken 
larger in order to reproduce experimental values of /^r, Am^ and S at the same time in 
the calculations based on the BS equation with the improved ladder approximation. 
(When we take tp = 0.04, the value of the running coupling in the infrared energy 
region becomes about a ~ 35.) 

Now let us look at what other quantities can be obtained from nyy^n^^. As for p 
meson and ai meson, we derive their masses and decay constants by fitting them to the 
calculated UyviQ^) — ^aa{Q^) using the pole saturated form of UyviQ'^) — ^aa{Q^)- 
Here, we use the following simplest version of the pole saturated form as a fitting 
function: 

nMQ^) - nMQ%, = - + + jg^- (5.i3) 

Resultant best fitted values are listed in table El together with the experimental val- 
ues. From the results shown in this table, we can see that the masses and the decay 
constants of p and ai mesons are not so much dependent on the infrared cutoff param- 
eter tp compared with tp dependence of Am^ and 5*. When we choose tp = 0.04, for 
which the values of Am^ and 5* are predicted in good agreement with experiments. 
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tp 


0.04 


0.05 


0.1 


0.2 


0.3 


0.5 


Exp. 


Mp (MeV) 


730 


710 


716 


645 


613 


612 


770 - 772 


/p (MeV) 


150 


146 


151 


150 


150 


151 


154 


Ma, (MeV) 


1186 


1146 


1012 


908 


861 


859 


1190 - 1270 


fa, (MeV) 


115 


109 


119 


118 


119 


123 


144 



Table 3: Resultant best fitted values of masses and decay constants of p and oi mesons 
for several values oitp- Experimental values for each quantities are also listed. Exper- 
imental values of Mp, fp and Ala^ are given in Ref. and fa, is given in Ref. |37| . 



we obtained 

Mp = 730 MeV, (5.14) 

fp = 150 MeV, (5.15) 

Ma, = 1186 MeV, (5.16) 

fa, = 115 MeV.. (5.17) 

These are also in good agreement with experiments. This confirms the rehabihty of 
calculations in the present analysis based on the improved ladder BS equation with 
small infrared cutoff parameter {tp ^ 0.04). 



6 Summary and Discussions 

We have calculated the 7r~^ — 7r° mass difference Ami same footing as the pion 

decay constant f^^ and the QCD S parameter through the vector and the axial-vector 
current correlators. For the calculations we solved inhomogeneous BS equations for 
the vector and axial-vector vertex functions together with SD equation for the quark 
mass function within the improved ladder approximation in the Landau gauge. We 
also obtained the masses and the decay constants of p and ai mesons by fitting them 
to UvviQ'^) — n^A(<5^) using the pole saturated form. 

We showed that all of these quantities are simultaneously fit in agreement with 
the experiments when we take the IR cutoff parameter tp = 0.04, which corresponds 
to taking 

AQCD^ 724MeV. (6.1) 

This parameter choice of tp is fairly smaller than the ones investigated in the previous 
works 1201 1211 122 12^1 , which corresponds to 

AQCD^ 500MeV. (6.2) 

Either value of Aqcd above is substantially larger than that in the MS scheme |24j : 

Ag),j3 = 300 - 450 MeV . (6.3) 

Our results imply that the running coupling of QCD in the IR region should be taken 
larger than that considered so far, when we calculate not only /^r but also the QCD 5* 
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parameter and Am^ in the framework of the BS and the SD equations in the improved 
ladder approximation. Apparently, choosing the IR cutoff parameter corresponds to 
exploiting the "scale ambiguity" jSHj of Aqcd in the SD and BS equations as was 
recently emphasized in Ref. jSEI which proposed an interesting method to resolve the 
scale ambiguity in the the calculation of fj^. 

In order to establish the approach of the SD and the BS equations in the improved 
ladder approximation, it is certainly desirable to apply the method of Ref. [261 
case, namely the calculations of Am^ and QCD S parameter as well as /^r. This will 
be done in future work. 

The success of the analysis based on the BS and the SD equations in the real-life 
QCD in this paper well motivate us to apply this method for studies of other strong 
coupling gauge theories. One of the most interesting examples of such is the large Nf 
QCD (see, e.g., Refs. 1^11^311011111113113111), which is a QCD with large number 
of massless flavors (but not too large number as to destroy the asymptotic freedom). 
It is expected that the chiral symmetry gets restored for a certain large number of 
massless flavors, which was in fact confirmed by the lattice simulations jH]. 

It was argued HH Ej based on the effective field theory of HLS [ZllHl Ej that 
this chiral restoration of the large Nf QCD is accompanied by the massless vector 
meson degenerate with the NG boson (pseudoscalar meson) as the chiral partners 
("Vector Manifestation" of chiral symmetry). It was further argued |H| that the HLS 
model even in the real-life QCD behaves as a little Higgs model |2j with two sites 
and two links, whose locality of theory space forbids the quadratic divergence in Am^ 
which is nothing but the (mass)^ of the Higgs in the framework of the Little Higgs 
models. In the large Nf QCD Am^ becomes small compared with near the chiral 
restoration point, which may suggest that the large Nf QCD may be considered as a 
UV completion of the Little Higgs model. 

On the other hand, we found in the previous paper ^Hj , by the explicit calculation 
based on the homogeneous BS and the SD equations in the improved ladder approx- 
imation, that masses of the scalar, vector, and axial-vector mesons in the large Nf 
QCD are proportional to /^r and vanish at the chiral restoration point, which implies 
somewhat different manifestation than the Vector Manifestation in the HLS model. 
(This behavior was also conjectured in Ref. jlHI-) 

Thus it would be very useful to clarify the situation to calculate the S parameter 
and Am^ in large Nf QCD which have not yet been calculated directly from QCD. 
In the forthcoming paper we shall calculate Am^ and S parameter in large Nf QCD 
by the same method presented in this paper. 
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